Interplay between Josephson effect and magnetic interactions in double quantum dots 
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We analyze the magnetic and transport properties of a double quantum dot coupled to supercon- 
ducting leads. In addition to the possible phase transition to a 7r state, already present in the single 
dot case, this system exhibits a richer magnetic behavior due to the competition between Kondo and 
inter-dot antiferromagnetic coupling. We obtain results for the Josephson current which may help 
to understand recent experiments on superconductor-metallofullerene dimer junctions. We show 
that in such a system the Josephson effect can be used to control its magnetic configuration. 

PACS numbers: 74.50.+r,73.63.-b,75.20.Hr,73.21.La 



-a 
c 

o 
o 



> 

in 
m 
m 
O 

o 

-i— > 



i 

C 

o 
o 



X 



Quantum dot (QD) devices provide a unique opportu- 
nity to study the interplay between different basic quan- 
tum phenomena. Thus, for instance, great advances in 
the understanding of Kondo physics have been achieved 
since the observation of the Kondo effect in semiconduct- 
ing quantum dots [l| . More recently double quantum dot 
(DQD) structures have been proposed for studying the 
competition between the Kondo effect and the inter-dot 
antiferromagnetic coupling |2|. An additional interest- 
ing ingredient is introduced when these systems are con- 
nected to superconducting electrodes In this case 
the electron pairing in the leads appears as a competing 
mechanism to both the Kondo and other type of mag- 
netic interactions that could be present. For a single 
quantum dot placed between two superconductors this 
competition can lead to a suppression of the Kondo ef- 
fect and the appearance of an unscreened magnetic mo- 
ment, corresponding to a quantum phase transition to 
the so-called 7r-state with a reversal of the sign of the 
Josephson current 0, 0. On the experimental side, 
great progress in the physical realization of these sys- 
tems is being achieved by structures consisting of nan- 
otubes or fullerene molecules attached to metallic elec- 
trodes |E 0- In R°f- electron transport through 
superconductor- mctallofullerene molecules (Gd@Cg2)- 
superconductor junctions was analyzed. Strong features 
associated with superconductivity were observed for the 
case of junctions containing a molecular dimer. As 
pointed out in |7|, the observed non-monotonic depen- 
dence of the low bias current as a function of temperature 
could be related to a change in the magnetic configura- 
tion of the Gd atoms. DQD systems coupled to super- 
conducting electrodes have been theoretically analyzed 
in Refs. However, these works considered geome- 

tries and ranges of parameters which do not correspond 
directly to the situation in the experiments mentioned 
above. 

In this Letter we provide an analysis, based on exact 
diagonalizations and mean field slave boson techniques, 
of the interplay between the Josephson effect, Kondo cor- 
relations and antiferromagnetic coupling in S-DQD-S sys- 
tems. Like in the single S-QD-S case we identify phases 
in which the sign of the Josephson coupling is reversed. 



The situation in the S-DQD-S system is however richer 
from the point of view of its magnetic configuration. We 
show that when the system is coupled to localized spins 
as in the experimental situation of Ref. their relative 
orientation can be influenced by the Josephson current 
through the device. We claim that these properties pro- 
vide a way to control the magnetic configuration of such 
a nanoscale system. 

The system depicted in Fig. 1 consists of two coupled 
quantum dots in series, placed between two supercon- 
ducting electrodes. The electronic degrees of freedom 
are represented by a double Anderson model with a sin- 
gle spin-degenerate level on each QD. The corresponding 
Hamiltonian is given by 

H el = H L + H R + e ia h ia + U ^ n^n il 

+H 12 + H XL + H 2R , (1) 

where the index i = 1,2 identify each QD; the terms 
Hr, and H R describe the uncoupled leads as BCS super- 
conductors characterized by a complex order parameter 
Ae"* L ' fl . H12 is the coupling term between the dots given 
by H\2 — J2cr^ 12 ^ia^a + h.c. The last two terms cor- 
respond to the coupling between the dots and the elec- 
trodes, i?i£(2H) = Y,k a t L,Rc\ {2)!y c kL ( R - )a +h.c. The 
Coulomb interaction within each dot is described by the 
U term. The coupling of the QD's with two localized 
magnetic moments, which in the experiments of Ref. 
are provided by the Gd impurities, can be modeled by 
an additional term in the Hamiltonian of the form 



Hint = JS\.a\ + JS2-&2 



(2) 



where a\ t 2 are the electronic spin operators in the dots, 
while S\^2 denote the localized spins. We assume, in ac- 
cordance to Refs. that the magnetic coupling, 
J, is much weaker than the other energies involved in 
the problem, which allows to introduce its effect as a 
perturbation in a second stage. There exists some con- 
troversy concerning the sign of J in Gd@C82 although 
recent studies suggest that this coupling is antiferromag- 
netic [lj}. The actual sign of J is, however, not essential 
for the main effects discussed below. These studies also 
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FIG. 1: (color online) Schematic view of the studied structure 
(a) and involved energies (b). This system should model two 
fullerenes doped with Gd atoms in contact with two super- 
conducting reservoirs (L, R). 



suggest a large magnetic moment associated to the Gd 
impurities, which allows to consider Si t 2 m our model as 
classical. We will focus in the case of strong coupling 
between the QD's by taking £12 /A = 10, which roughly 
corresponds to the estimates [ll| for the experimental 
situation of Ref. 0. The results of this reference also 
suggest good coupling to the leads (i.e. tL,R > A). 

A first insight into this problem can be provided by 
analyzing the H e i ground state properties as a function 
of the model parameters. For this purpose we rely on an 
approximation consisting in taking the zero bandwidth 
limit (ZBWL) for the superconducting electrodes. The 
validity of this approach has been discussed for other su- 
perconducting junctions in Refs. [5|, Il4| . In this limit 
the Hilbert space of the S-DQD-S system is restricted 
to 4 4 states and H e i can be diagonalized exactly. In 
the superconducting case we distinguish four different 
ground states: the pure and tt states for which the 
energy as a function of the superconducting phase dif- 
ference <j) — 4>l — 4>r has a minimum at (f> = and tt 
respectively; and two intermediate phases, which are de- 
signed as 0' and tt' depending of the relative stability 
of each minima Fig. 2 illustrates the (e, U) phase 
diagram for two different values of t R = t L : 2 A (panel 
(a)) and 2.5A (panel (b)). We show only the range of e 
which corresponds to a charge per dot varying between 
and 1, where the transition to the tt state can take 
place 12]. As shown in Fig. 2, for t L = t R = 2A all 
phases 0, 0', tt' and tt appear at the transition region. A 
more detailed understanding of the ground state proper- 
ties is provided by analyzing non-local spin correlation 
functions of the form < a^a^ >. We choose the line in 
the phase diagram which corresponds to a large intradot 
Coulomb interaction, U — 800A, and show the evolution 
of these correlation functions with e in Fig. 2 (c) together 



with the occupation numbers nj^ for each dot. The ap- 
pearance of a 1/2 magnetic moment for the full S-DQD-S 
system is signaled by the broken symmetry ^ n j . The 
function < a\di> measures the correlation between the 
electron spins in the two dots. As can be seen in the 
middle panel of Fig. 2 (c) it evolves continuously from 
to —3/4, the latter value corresponding to a complete 
antiferromagnetic (AF) correlation. It is worth noticing 
that this AF tendency is more pronounced in the super- 
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FIG. 2: (color online) (U, e)-phase diagram for ti2 = 10A, 
tL — tu — 2A (a) and tL = tn = 2.5A (b) indicating the 0, 
0', 7r' and tt regions. On (c) we show the e dependence along 
the line U = 800 A indicated in (a) of (from top to bottom): 
the occupation number for spin up and spin down electrons; 
the interdot spin correlation function, and between the DQD 
spin and the electrodes. The dashed (solid) lines corresponds 
to the normal (superconducting) state. 



conducting case as the presence of the superconducting 
gap reduces the number of low energy excitations capable 
to screen the spin in the dot region, i.e. it leads to a par- 
tial suppression of the Kondo correlations. In the present 
range of parameters, with a strong interdot hopping, the 
Kondo regime corresponds to the formation of a spin sin- 
glet between an electron in the bonding state of the DQD 
"molecule" and the electrons in the leads. This correla- 
tion is reflected in the behavior of < &l(R)&i(2) >j which 
becomes increasingly negative in the Kondo regime. As 
can be observed in the lower panel of Fig. 2 (c), Kondo 
correlations are strongly suppressed by superconductiv- 
ity compared to the normal case. In fact, it is in the range 
of parameters corresponding to the Kondo region in the 
normal state where the tt state appears. The normal 
state itself exhibits a transition from the Kondo to the 
AF regime for — e ;§> ti2, signaled by < <7l(r)<?i(2) >^ 
and < f?i<?2 >^ —3/4. This transition roughly corre- 
sponds to the situation where the AF coupling between 
the dots ~ ti 2 /U becomes larger than the AF coupling of 
the dots with the leads ~ 2t\ R /\e — t\%\. If the hopping 
to the leads is increased, the suppression of the Kondo 
effect by superconductivity becomes less effective. As 
shown in Fig. 2 (b), for t L = t/R = 2.5A the system ex- 
hibits only three phases (0, 0' and tt') within the range 
of parameters considered. Further increase of tL,R would 
lead to a complete suppression of the tt' and 0' phases. 

In order to go beyond the ZBWL approximation and 
include the finite bandwidth of the electrodes we use an 
appropriate slave-boson representation of Hamiltonian 
Q. In Ref. the U — > oo mean field slave-boson 
approach was used to study the single QD system 
with superconducting electrodes. However, in order to 
describe the main features of these systems including the 
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possibility of unscreened magnetic moments it is neces- 
sary to use the more general representation of Ref. 0] , 
which is valid for finite values of U and allows for pos- 
sible magnetic solutions ^t|- Following Ref. ^(| the 
auxiliary Bose fields are designed by ii (empty state), 
Pi a (single occupied state corresponding to spin a) and 



di (double occupied state) and we define the operator 
z la = {l-dl-f^^ie^+P^il-el-f^)- 1 / 2 , 
where i = 1,2 denotes the two different QDs. In the 
enlarged space the Hamiltonian has the form 



t L(R){z\^ 2 - }a f\ i 2- )cr C kL ( R - )cJ + h.C.) 

icr i a k,a 

- a i(^A + ^di +^p\ a Pia ~ 1) - Pivifiaf™ ~ v\ a Vta ~ <$ <k) (3) 
i tr ia 

I 



where the fi a are fcrmionic operators and oij, f3i a are 
the Lagrange multipliers corresponding to the constraints 
e\e t + dUi + J2ctPLpw = 1 and j\J ia = p\ a p la + dftdi. 
The particular definition of the z\ a operators warrants 
that the exact solution in the U — > limit is recovered 
|16| |. Within a mean field approximation we replace the 
Bose operators in Eq. @ by their mean values e^, di, 
Pic and Zi a . In this approximation the Hamiltonian pa- 
rameters are renormalized according to l a i — ej — @i a , 
h2a = t 12 z la z 2 o and i L{R)rT = t L ( R) z 1{2)a . The mean 
values of the Bose operators must be determined self- 
consistently by minimizing the effective action |16| . 

Both the current through the dots I = ie/h^2 a t\2 a [< 
flaha > - < flaha >] and their occupations n ia =< 
r a ha > can be calculated using standard Green func- 
tion techniques ^^|. We have solved numerically the 
mean field equations and computed the Josephson cur- 
rent through the DQD system for a certain set of param- 
eters. In Fig. 3 we show the current-phase relation for 
U = 800A, ti2 = 10A and two different values of the pa- 
rameter T L , R = Trt 2 L R p LtR (E F ), where Pl,r( e f) is the 
electrodes normal density of states at the Fermi energy. 
We see that for Tl = T R = 2. 25 A the system evolves 
from the to the ir state as e varies between ~ — 5A and 
~ — 2A going back to the state for e ~ 7A. In the case 
of a larger coupling, Tl,r = 4A, the pure ir state is never 
reached, in good qualitatively agreement with the results 
obtained within the ZBWL. It is worth noticing that the 
occurrence of the ir state r equires A to be larger than 
an energy scale ~ y/TU/2 exp (— 7r|e — ti 2 |/2r), associ- 
ated with the Kondo effect of the singly occupied bond- 
ing level, which for ti 2 /T > 1 can be much smaller than 
the effective Kondo temperature estimated for the nor- 
mal DQD system |l9|. 

Let us now analyze how the behavior of the electronic 
system could influence the configuration of the localized 
spins by means of their magnetic coupling given by Eq. 
(J2J) . One would expect that the appearance the magnetic 
7r-state could give rise to a change in this configuration 
with respect to the case of normal electrodes. In fact, for 




FIG. 3: (color online) The current-phase relation for t\2 
10A, U = 800A, F L = F R = 2.25A (upper panel) and F L 
Fr — 4A (lower panel) 



small J the total energy can be expanded as 
E(h 1 ,h 2 )^E(Q,0)+ J2 O/Ai + S 



(4) 



where = JS^^, S^. z being the z-component of the 
localized spin. Notice that in general, for non-magnetic 
situations, only the quadratic correction appears as in 
the well known RKKY interaction [2(]]. For the range 
of parameters where the 7r states appear this correction 
would be positive leading to a AF configuration of the lo- 
calized spins. However, the broken symmetry, nj ^ n\, 
in the 7r state gives rise to non-vanishing linear correc- 
tions which favor the parallel (F) configuration. This is 
illustrated by the insets in the upper panel of Fig. 4 
which show the behavior of the total energy as a func- 
tion of h = \hi\ = at <f> — and <f> — tt respectively 
in the 0' region of the phase diagram. A simple image 
of this effect is that the unscreened magnetic moment 
appearing in the 7r-state acts as a local magnetic field 
which tends to align the localized spins. The full phase 
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FIG. 4: (color online) Energy (upper panel) and current- 
phase relation (lower panel) corresponding to the parallel (F, 
dashed line) and antiparallel (AF, full line) configuration of 
the localized spins. We have chosen U = 800A, ti2 = 10A, 
\hi\ = \h 2 \ = 0.25A and T L = F R = 2.25A. The insets show 
the behavior of the energy as a function of \hi,2 1 for the F and 
AF configuration at (j> — and </> = n. 



dependence of the total energy for finite J in the F and 
the AF configurations is depicted in the upper panel of 
Fig. 4. As can be observed the range of stability of the 
7r state is increased in the F configuration, which in turn 
has a noticeable effect in the Josephson current (lower 
panel of Fig. 4). Thus, our results predict that the sys- 
tem switches from the AF to the F configuration as the 



superconducting phase is sweeped from to tt. They 
also suggest that in the transition region the F configu- 
ration is a metastable solution which could give rise to an 
hysteretic behavior as a function of the superconducting 
phase difference. Although a direct comparison with the 
data of Ref. ,?] is not possible since the experiment was 
performed under non-equilibrium conditions, our results 
tend to support that the non-monotonic behavior of the 
low bias current as a function of temperature may be in- 
deed related to a change of the magnetic configuration 
of the Gd atoms. A more direct test of our predictions 
would require the measurement of the supercurrent in a 
phase biased situation. 

In conclusion we have studied the magnetic and super- 
conducting properties of a S-DQD-S system. We have 
shown that it can exhibit a quantum phase transition to 
a 7r state with an unscreened 1/2 magnetic moment in 
the dots region. When the system is coupled to local- 
ized spins as in the experimental situation of Ref. pj a 
transition from an AF to a F configuration can be in- 
duced by tunning the superconducting phase difference. 
These properties illustrate the possibility of controlling 
the magnetic configuration at the nanoscale by means of 
the Josephson effect. 
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